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Abstract 

Higgs boson production via weak boson fusion is sensitive to the tensor structure of the HVV 
iy = W, Z) couplings, which distinguishes loop induced vertices from SM expectations. At the 
CERN Large Hadron Collider this information shows up most clearly in the azimuthal angle cor- 
relations of the two forward and backward quark jets which are typical for weak boson fusion. We 
calculate the next-to-leading order QCD corrections to this process, in the presence of anomalous 
HVV couplings. Gluon emission does not significantly change the azimuthal jet correlations. 



1 



Introduction. The production of Higgs bosons in the weak boson fusion (WBF) process 
will provide a direct and highly sensitive probe of HWW and HZZ couplings at the CERN 
Large Hadron Collider (LHC) 

mm 

. The determination both of the strength and of 
the tensor structure of these couplings is crucial for the identification of the produced boson 
as a remnant of the spontaneous symmetry breaking process which is responsible for W and 
Z mass generation. 

Within spontaneously broken, renormalizable gauge theories like the standard model 
(SM), this coupling originates from the kinetic energy term, (Z)^$)^(D'^$), of a scalar 
Higgs field, $, whose neutral component obtains a vacuum expectation value (vev), 
$° ^ {v + H)/\^. This replacement then leads to a characteristic coupling in the in- 
teraction Lagrangian, of the form HV^jV^ {V = W, Z). The existence of the vev is necessary 
to produce a trilinear HVV coupling at tree level: with f = all couplings to the gauge 
fields V contain two scalar fields, i.e., only HHV and HHVV couplings would be gener- 
ated. A trilinear HVV coupling may also be loop-induced, however. The SM H'j'y and Hgg 
effective couplings are an example: they are induced by W^-boson and/or top quark loops. 
Gauge invariance dictates a different tensor structure of these loop-induced couplings: the 
corresponding effective Lagrangian contains the square of the field strength, i.e. the lowest 
order loop-induced terms are of the form HV^^V^^ or HV^yV^^ , where V^'' = \e'^^P"Vpa 
denotes the dual field strength of the gauge field. 

The task of future Higgs experiments is, then, twofold: (i) to measure the overall strength 
of the HVV coupling, and (ii) to identify its tensor structure. One would expect a loop- 
induced coupling to be much smaller than the expected SM HVV coupling strength. How- 
ever, the measurement of WBF rates alone will not be sufficient to establish H as being 
related to spontaneous symmetry breaking: to give just two examples, the loop-induced 
couplings might be substantially enhanced by additional non-SM particles in the loop or by 
the existence of multiplets of large weak isospin which couple strongly to H . Or a particular 
LHC signature may be strongly enhanced by a much larger H decay branching ratio than 
in the SM. A confirmation that the HVV coupling has tree level strength is, thus, ambigu- 
ous: a clear identification of the Higgs boson also requires the identification of the tensor 
structure of the HVV vertex. 

It was pointed out some time ago that the azimuthal angle correlations of the two quark 
jets in the weak boson fusion process qQ qQH provide tell-tale signatures for the tensor 
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structure of the HVV couplings : the SM expectation is for a flat distribution, while the 
loop-induced couplings lead to a pronounced dip at azimuthal separations (pjj of the two 
tagging jets of 90 degrees for a HV^yV^" coupling and at and 180 degrees for the CP 
violating HV^^V^^ vertex. Observation of the tagging jets is crucial for isolating the WBF 
process from backgrounds and, therefore, their distributions will be available for all WBF 
samples. Also, signal to background ratios for WBF processes are expected to be very good 
within the SM, exceeding the 1:1 level for wide ranges of the Higgs boson mass 

The analysis of Ref. [6| was performed at leading order (LO) in QCD. This means that 
additional gluon emission, which might lead to a de-correlation of the tagging jets, was 
ignored in the analysis. Subsequently it was argued ~\ that such de-correlation effects play 
an important role in a related process, gg — > Hgg, when the two tagging jets are widely 
separated in rapidity, which is a typical requirement for WBF studies. In this Letter we 
analyze this question, by calculating the tagging jet distributions in next-to-leading order 
(NLO) QCD, for the production of a scalar H via WBF with an arbitrary tensor structure 
of the HVV vertex. If de-correlation is important, it should show up in the form of large 
radiative corrections at NLO. We use the term "Higgs boson" as a generic name for the 
produced scalar in the following. 

The NLO calculation. Our calculation is an extension of the NLO QCD corrections 



BJaQ 



For the 



for the SM WBF processes qQ qQH (and crossing related ones) 
total cross section these corrections have been known for over a decade Recently, we 
have recalculated them by developing a NLO parton level Monte Carlo program which 
provides the flexibility to calculate arbitrary distributions at NLO, such as the azimuthal 
angle correlations that we are interested in here. 

uses a SM vertex function, T^'^{qi, ^2) = —;^g^^ for the HVV 
vertex in Fig. ^ Here we need to generalize this vertex to the most general structure 
compatible with Lorentz invariance. Taking into account that the quark currents in Fig. Q 
and for the corresponding gluon emission processes are conserved, all terms proportional to 
qi or 52 niay be dropped, and the most general HVV vertex may be written as 

T^%qi, q2) = a,{qu q^) g^" + a^iqi, ^2) [gi ■ q2g''' - q'2qV\ + a,{q,, q^) e^^^^-q^.q^^ . (1) 

Here qi and q2 are the four- momenta of the two weak bosons, and the aj(gi, ^2) are Lorentz- 
invariant form factors, which might, for example, represent scalar loop integrals in a per- 
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FIG. 1: Feynman graphs contributing to qQ qQH at (a) tree level and (b) including virtual 
corrections to the upper quark line. The momentum labels and Lorentz indices for the internal 
weak bosons correspond to the vertex function of Eq. (^. 

turbative calculation. It is straightforward to implement the general vertex of Eq. into 
our NLO QCD Monte Carlo: the virtual amplitude of Fig. Q is proportional to the Born 
amplitude, A^Born, irrespective of the structure of the HVV vertex. Thus, all amplitudes 
reduce to a simple contraction of quark (or quark-gluon) currents with the vertex function 
of Eq. (^. These currents, and their contractions, are evaluated numerically, using the am- 
plitude formalism of Ref . [11| . All other aspects of the present NLO calculation are handled 
as in Ref. 0], except that we do not simulate any Higgs boson decays in the following. Fac- 
torization and renormalization scales are fixed to fip = fJ'R = Qi for QCD corrections to the 
first or second quark line in Fig. ^ Here Qi and Q2 are the virtualities of the exchanged 
weak bosons. We use CTEQ6M parton distributions Q with as{Mz) = 0.118 for all NLO 
results and CTEQ6L1 parton distributions for all leading order cross sections. 

Anomalous couplings and form factors. While the g'^'^-teim in the vertex function 
corresponds to a SM Higgs coupling, the anomalous coupling terms 02 and 03 can be re- 
lated to higher dimensional operators in an effective Lagrangian. They first appear at the 
dimension-5 level and may be written as 



HWW HWW 

^S^EZ^^Z^^ + ^-^HZ.^Z^^ , (2) 

2A5e ^ 2A50 ^ ^ ' 

where the subscript e or o refers to the CP even or odd nature of the individual operators. 
In our discussion we will neglect possible contributions from if77 and if 7Z couplings which 



can appear in SU{2) x f/(l) invariant formulations jisl . Q|. The precise mix of HWW, 
HZZ, HZ'-f and H'-f-f contributions is quite irrelevant for the observable azimuthal angle 
distributions, as long as we do not consider interference effects between SM and anomalous 
vertices, and it will not affect our conclusions about the size of NLO corrections. For 
simplicity we therefore set Oi = for the anomalous coupling case and choose relative 
contributions from WW and ZZ fusion as in the SM, by taking g^J^^ = g^^^ = 1, 
9be^^ = 9bo^^ = l/cos^6'vy, and by using either ~480 GeV, A50 = 00 for the CP even 
case or A50 ~480 GeV, Ase = 00 for the CP odd case, which roughly reproduces SM rates 
for a scalar mass of mn = 120 GeV. 

The effective Lagrangian of Eq. Q produces couplings 

a2(gi, q.) = --^^sT"^ , a3(g„ g,) = -^gg^^ (3) 

for the HWW vertex, and 

a2(gi, q2) = -^g"e^^ , as{q,, ga) = -^g^o^^ (4) 

for the HZZ vertex. In general, the Oj are form factors which are expected to be suppressed 



once the momentum transfer, \J —qf-, carried by the virtual gauge boson reaches the typical 
mass scale, M, of the new physics which is responsible for these anomalous couplings. Below 
we use the simple ansatz 

a,(gi,g2) = a.(0,0)^^-^^^-^ (5) 

for discussing the consequences of such form factor effects. 

Results: The typical signature of a weak boson fusion event at the LHC consists of the two 
quark jets (tagging jets) and the Higgs decay products. The tagging jets tend to be widely 
separated in rapidity, with one quite forward (typical pseudorapidity of 3 to 4) and the second 
one backward, but frequently still located in the central detector (pseudorapidity below 2.5). 
Various Higgs decay modes have been considered in the literature for WBF, H WW [l|, 
H ^ TT 2], and if 77 being the most promising ones. While optimized event selection 
varies, in particular for the decay products, the cuts on the tagging jets are fairly similar in 
all analyses. Since here we are interested in the QCD features of WBF events, which do not 
depend on the Higgs decay mode, we perform our NLO analysis without simulating Higgs 
decays, and we only impose typical WBF cuts on the tagging jets. 




FIG. 2: Normalized transverse momentum distribution of the hardest jet for the SM Higgs boson 
(soHd red hne) and a scalar H of mass niH = 120 GeV with CP even anomalous coupling 02(^1, ^2)- 
The dash-dotted curves correspond to different form factor scales M = 100, 200, 400 GeV in Eq. Q 
and 02 = const, (blue curves) at NLO. LO curves are shown by the dashed lines and differ very 
little from the NLO results. 

In order to reconstruct jets from the final-state partons, the fcy-algorithm as described 
in Ref. [l^ is used, with resolution parameter D = 0.8. In a given event, the tagging jets 
are then defined as the two jets with the highest transverse momentum, ptj, with 

PTj> 20 GeV, |%|<4.5. (6) 

Here yj denotes the rapidity of the (massive) jet momentum which is reconstructed as the 
four- vector sum of massless partons of pseudorapidity Ir^l < 5. Backgrounds to weak-boson 
fusion are significantly suppressed by requiring a large rapidity separation of the two tagging 
jets. This motivates the final cut 

^Vjj = Ivji - %2l > 4 , Vj, ■ < , (7) 

which includes the requirement that the two tagging jets reside in opposite detector hemi- 
spheres. 
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The structure of the HVV couphng affects the production dynamics of H and we can 
expect significant deviations in jet observables if, instead of the SM, anomalous couphngs 
describe the vertex of Eq. ([Q). One example is shown in Fig.|2l where transverse momentum 
distributions, da / dp^jimsix) , are compared between the SM (solid line) and the CP even 
coupling 02(^1,(72)5 with different form factor scales M in Eq. Here, ^^■^(max) is the 
maximum pt of the two tagging jets. Only the shape of the distribution is considered, 
since the rate can always be adjusted by multiplying the anomalous couplings by a constant 
factor. Also, we should note that a CP odd coupling leads to very similar curves for a given 
form factor scale. In all cases we show the LO expectations (dashed lines) together with the 
NLO results: QCD corrections are of order 10%, typically, and well under control. 

One finds that anomalous HVV couplings generally lead to harder pj- spectra of the 
two tagging jets. Since the anomalous Lagrangian in Eq. (j2)) couples the Higgs boson to 
weak boson field strengths, transverse polarizations of the incident VV pairs dominate the 
anomalous case, while longitudinal VV fusion is responsible for SM Higgs production. A 
telltale sign of transverse vector boson fusion is the more central and, hence, higher pt 
production of the tagging jets. This effect is enhanced by the momentum factors in the 
HVV anomalous vertices. 

While the changed transverse momentum distributions in Fig. |21 could be used to rule 
out the SM, the reverse is not readily possible: a jet transverse momentum distribution 
compatible with SM expectations might be faked by anomalous couplings and a judiciously 
chosen form factor behavior of the coefficient functions 02 or 03 in Eq. The different scale 
choices in Fig. El demonstrate this effect: a low form factor scale of M = 100 GeV or slightly 
lower would be difficult to distinguish from the SM expectation and one can certainly find 
a functional form of the form factors which reproduces the SM within experimental errors. 

A much better observable for distinguishing the different tensor structures of the HVV 
vertex is the azimuthal angle correlation of the two tagging jets, da/d(j)jj jj^. Here (pjj is the 
azimuthal angle between the two tagging jets. The corresponding distributions are shown 
in Fig. inifor the SM (solid line) and for the same choices of form factors as before. The dip 
at (f)jj = 90 degrees for the CP even coupling and the suppression at and 180 degrees for 
the CP odd coupling are clean signatures which only depend on the tensor structure of the 
couplings and not on the precise dynamics which is responsible for the form factors. The 
remaining form factor dependence is very small and can be explained by kinematic effects 
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FIG. 3: Normalized azimuthal angle distribution, 1/cj da/d(f)jj where (f)jj is the azimuthal angle 
separation of the two tagging jets. NLO (solid and dot-dashed) and LO results (dashed lines) 
are shown for niH = 120 GeV in the SM (red curves) and (a) for a CP even anomalous coupling 
'^2('i'i, 92)5 (b) for a CP odd anomalous coupling a^{qi-,q2) with form factor scales M = 100, 200, 
400 GeV and (blue curves) M = 00. 

related to the higher average jet transverse momentum for big form factor scales, M: at 
small (f)jj two high pt jets recoil against the H scalar, resulting in an increased invariant 
mass of the event compared to the situation with two back-to-back jets. This leads to a 
more asymmetric (f)jj distribution for high form factor scales. 

The pronounced dip at 90 degrees, which is characteristic of the CP even coupling, is 



also found in Hjj production via gluon fusion [17|, at LO. This is not surprising because, 
in the large top mass limit, the Hgg vertex can be described by an effective Lagrangian 
proportional to H G'^^^G"''^'^ , which exhibits the same field strength squared behavior and 
hence the same tensor structure as the CP even HVV coupling in Eqs. ()1I2|) . Since the two 
tagging jets are far apart from each other, separated by a large rapidity gap of 4 units of 
rapidity or more, this LO behavior may be significantly reduced by gluon radiation when 
higher order QCD corrections are taken into account. Such de-correlation effects have been 

n Q 

studied for dijet events at the Tevatron [18]. For Hjj production via gluon fusion, Odagiri [7|] 
has argued that the dip structure is largely washed out by additional gluon emission between 
the two tagging jets. 
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FIG. 4: Higgs mass dependence of the azimuthal angle separation of the two tagging jets. In 
(a) the normaUzed azimuthal angle distributions are shown at LO (dashed lines) and NLO (solid 
lines) for Higgs masses of mu = 120, 200, 500 GeV and a constant CP even anomalous coupling 
02. Corresponding K-factors are shown in (b). 

Our NLO calculations show that such de- correlation effects are irrelevant for weak boson 
fusion, where t-channel color singlet exchange severely suppresses gluon radiation in the 
central region. The LO and the NLO curves in Fig. IHl are virtually indistinguishable. In 
order to better exhibit the size of NLO QCD effects for the WBF case, we show, in Fig. IHa) 
the azimuthal angle correlations for a pure CP even anomalous coupling for three different 
Higgs masses, vnu = 120, 200 and 500 GeV. Only small changes are visible when going from 
LO (dashed lines) to NLO (solid lines). The differences between LO and NLO are smaller 
than kinematical effects that can be induced by cuts on the Higgs decay products or by 
variations of the Higgs boson mass. 

The small to modest size of the QCD corrections is quantified in Fig. El^b) where the K 
factor for the distribution is shown, which is defined as 



The i^-factor is below 1.4 even in the dip region, where the cross section is severely 
suppressed. Virtually identical results hold for the CP-odd case. Clearly, the characteristic 
azimuthal angle distributions of the jets in WBF are not affected in any significant way by 
NLO QCD corrections. 



Conclusions: We have performed a first calculation of the NLO QCD corrections to 
Higgs boson production via WBF in the presence of arbitrary anomalous HVV {V = W, Z) 
couplings. Anomalous couplings lead to characteristic changes in the azimuthal angle corre- 
lation of the two tagging jets in weak boson fusion events at the LHC, which provides for a 
very sensitive test of the tensor structure of the HVV couplings of the Higgs boson or of any 
other scalar with sufficiently large production cross section in WBF [Q]. We have shown by 
explicit calculation that these azimuthal correlations are not washed out by gluon emission, 
at NLO QCD, even though the tagging jets are widely separated in rapidity. This behav- 
ior can be understood as a consequence of t-channel color singlet exchange in WBF which 
severely suppresses the central gluon radiation which might cause tagging jet de-correlation. 

Acknowledgments 

This research was supported in part by the University of Wisconsin Research Commit- 
tee with funds granted by the Wisconsin Alumni Research Foundation and in part by the 
U.S. Department of Energy under Contract No. DE-FG02-95ER40896. 



[1] D. Rainwater and D. Zeppenfeld, Phys. Rev. D 60, 113004 (1999) [Erratum-ibid. D 61, 
099901 (2000)] iarXiv:hep-p h/9906218 ; N. Kauer, T. Plehn, D. Rainwater and D. Zeppenfeld, 
Phys. Lett. B 503, 113 (2001) arXiv :hep-ph/0012351 ; C. M. Buttar, R. S. Harper and 
K. Jakobs, ATL-PHYS-2002-033; K. Cranmer et al., ATL-PHYS-2003-002 and ATL-PHYS- 
2003-007; S. Asai et aL, ATL-PHYS-2003-005. 

[2] D. Rainwater, D. Zeppenfeld and K. Hagiwara, Phys. Rev. D 59, 014037 (1999) 
|arXiv:hep-p h/9808468 ; T. Plehn, D. Rainwater and D. Zeppenfeld, Phys. Rev. D61, 093005 
(2000) [arxiV:hep-ph/9911385| ; S. Asai et al., ATL-PHYS-2003-005. 

[3] D. Rainwater and D. Zeppenfeld, JHEP 9712, 005 (1997) |arXiv:hep-ph/9712271| ; K. Cran- 



mer, B. Mellado, W. Quayle and S. L. Wu, |arXiv:hep-ph/0401088 
[4] S. Asai et aL, |arXiv:hep- ph/0402254l 

[5] D. Zeppenfeld, R. Kinnunen, A. Nikitenko and E. Richter-Was, Phys. Rev. D62, 013009 
(2000) |arXiv:hep-ph/0002036 ; D. Zeppenfeld, in Proc. of the APS/DPF/DPB Summer 



10 



Study on the Future of Particle Physics (Snowmass 2001) ed. N. Graf, eConf C010630, 
P123 (2001) |arXiv: hep-ph/0203123 ; A. Belyaev and L. Reina, JHEP 0208, 041 (2002) 
[arXiv:hep-ph/0205270, . 

[6] T. Plehn, D. Rainwater and D. Zeppenfeld, Phys. Rev. Lett. 88, 051801 (2002) 

|arXiv:hep-ph/0105325] . 
[7] K. Odagiri, JHEP 0303, 009 (2003) |arXiv:hep-ph/0212215| . 

[8] T. Han, G. Valencia and S. Willenbrock, Phys. Rev. Lett. 69, 3274 (1992) 
rarXiv:hep-ph/9206246 . 

[9] T. Figy, C. Oleari and D. Zeppenfeld, Phys. Rev. D 68, 073005 (2003) 1arXiv:hep-ph/0306109| . 
[10] E. L. Berger and J. Campbell, |arXiv:hep-ph/0403194| 

[11] K. Hagiwara and D. Zeppenfeld, Nucl. Phys. B 274, 1 (1986); K. Hagiwara and D. Zeppenfeld, 

Nucl. Phys. B 313, 560 (1989). 
[12] J. Pumplin et al., JHEP 0207, 012 (2002) |arXiv:hep-ph/0201195| . 
[13] W. Buchmiiller and D. Wyler, Nucl. Phys. B 268, 621 (1986). 

[14] K. Hagiwara, S. Ishihara, R. Szalapski and D. Zeppenfeld, Phys. Rev. D 48, 2182 (1993). 
[15] S. Catani, Yu. L. Dokshitzer and B. R. Webber, Phys. Lett. B 285 291 (1992); S. Catani, 

Yu. L. Dokshitzer, M. H. Seymour and B. R. Webber, Nucl. Phys. B 406 187 (1993); S. D. Ellis 

and D. E. Soper, Phys. Rev. D48 3160 (1993). 
[16] G. C. Blazey et al, |arXiv:hep-ex/0005012"| 

[17] V. Del Duca, W. Kilgore, C. Oleari, C. Schmidt and D. Zeppenfeld, Phys. Rev. 
Lett. 87, 122001 (2001) |arXiv:hep-ph/0105129| ; Nucl. Phys. B 616, 367 (2001) 
farXiv:hep-p h/0108030 . 

[18] A. H. Mueller and H. Navelet, Nucl. Phys. B 282, 727 (1987); V. Del Duca and C. R. Schmidt, 
Phys. Rev. D 49, 4510 (1994) |arXiv:hep-ph/9311290 ; W. J. Stirling, Nucl. Phys. B 423, 56 



(1994) (arXiv:hep-ph /9401266 ; S. Abachi et al. [DO Collaboration], Phys. Rev. Lett. 77, 595 
(1996) |arXiv:hep-ex/9603010| ; L. H. Orr and W. J. Stirling, Phys. Rev. D 56, 5875 (1997) 
|arXiv:hep-ph /9706529 ; J. Kwiecinski, A. D. Martin, L. Motyka and J. Outhwaite, Phys. 
Lett. B 514, 355 (2001) |arXiv:hep-ph/0105039| . 
[19] The dimension 5 language is appropriate for, e.g., an isosinglet scalar resonance H. For a 
Higgs doublet $ with a vev, the leading operators appear at dimension 6 level [l^ . 1^ and 
the couplings in Eq. are suppressed by an additional factor g^^^ ~ w/A. 



11 



